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I. INTRODUCTION 


Graphene, a monoatomic layer material made of carbon 
atoms arranged in a honeycomb lattice, is a true 2D crystal 
structure. It had been known for decades that graphite, the 
material used in pencil tips, is made of stacked layers of 
graphene held together by weak Van der Waals forces. 
Moreover, the electronic properties of graphene and its 
relativistic band structure have been the object of intensive 
theoretical study in connection with the graphite used in 
nuclear reactors, [1]-[4]. However, before its first 
experimental realization in 2004 by a group at the University 
of Manchester in 2004 [2], it was believed that such a 2D 
crystal would be unstable due to thermal and other crystal 
fluctuations and could, therefore, not be practically 
obtained [3]. Following its aforementioned experimental 
validation and due to its unique electronic properties, 
graphene has triggered huge research efforts on its various 
features including electronic, electromagnetic, mechanical 
and thermal properties. 


More recently, graphene plasmonics has become a 
research area of particular interest in electromagnetics, due 
to several unique properties [5]. Our group has developed 
applications involving nonreciprocity, which can be tailored 
towards nonreciprocal plasmonic devices such as plasmonic 
isolators [6]-[15]. One of the methods for generating 
nonreciprocity in graphene is through a static magnetic field. 
The resulting magnetoplasmons and their nonreciprocity are 
the focus of this paper, which aims at presenting the relevant 
material in a self-consistent and pedagogical fashion. A brief 
introduction on the fundamentals of graphene, its band 
structure and its conductivity model is presented in Sections 
II and II. Section IV includes introductory material on 


surface plasmons and magnetoplasmons in graphene, and 
finally, Section V focuses on nonreciprocal properties of 
graphene magnetoplasmons and their applications. 


II. FUNDAMENTALS OF GRAPHENE 


First, a brief overview of graphene carbon bonding and 
crystallography is presented. These concepts are then applied 
in the tight binding model to derive the energy band structure 
of graphene. Finally, the concept of graphene’s effective 
Dirac Hamiltonian is introduced, as a tool to determine the 
quantized energy levels in magnetically biased graphene. 

A. Carbon bonds in graphene 


The carbon atom, the sole constituent of graphene, has 
four electrons in its outer shell, called the valence electrons, 


occupying the 2s° and 2 p atomic orbitals. In the 
crystalline phase, these weakly-bound electrons give rise to 
2s,2p,,2p,and2p, orbitals, which are important for 
forming covalent bonds in carbon materials [16]. Since the 
energy difference between 25 and 2p orbitals in carbon is 


small compared to the energy of the chemical bonds, the 
electronic wave function for these four electrons can readily 
mix with each other so as to enhance the binding energy of 
the C atom with its neighboring atoms [16]. This mixing of 


2s and 2p atomic orbitals is called hybridization, whereas 
the mixing of a single 2s electron with n=1,2,3 2p 
orbitals is called sp” hybridization. 


In the graphene 2D crystal, the outer shell electrons of the 
carbon atom form sp” hybridizations. The 2s and 2p, 


and 2p, orbitals mix to form three sp’ orbitals in the xy 
plane (plane of graphene) with angles of 120 degrees, 
leaving behind an unaltered 2p, orbital normal to the xy 


plane, shown in Fig. 1. 
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Fig. 1. Graphical representation of the SP orbitals in graphene. 


In graphene’s honeycomb structure, these three sp? 
orbitals form strong longitudinal in-plane covalent bonds 
with the sp? orbitals of the neighboring C atoms, called O 


bonds as shown in Fig. 2. O bonds are responsible for 
rigidity and the robust mechanical properties of graphene. 


The remaining 2p, orbital of each C atom forms lateral 
covalent bonds with the 2p, orbitals of the neighboring 


orbital of the 


neighboring C atoms have a small overlap, 7 bonds are 
weaker compared to the O bonds. 7 bonds are responsible 
for the electronic properties of graphene. 


atoms called 7 bonds. Since the 2p, 
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Fig. 2. Formation of ©O and 7 bonds between neighboring carbon 
atoms in graphene. 


B. Graphene crystallography 


Graphene’s honeycomb lattice is shown in Fig. 3 with 
solid lines showing the bonds between neighboring C atoms. 


The vectors 
= 5 (3% + V37) and a, = 5 (38 — V37) 


are the direct lattice vectors, where a % 1.42 A isthe inter- 
atomic distance. Note that the honeycomb structure is 
composed of two staggered triangular sub-lattices shown 


2 


with different colors in Fig. 3. The structure may also be 
seen as a triangular lattice with a 2-atom basis, represented 
as A and B in Fig. 3. Each unit cell of the honeycomb lattice 
contains 2 C atoms (6 carbon atoms at the corners of each 
hexagon, each atom shared by 3 hexagonal unit cells). 





Brillouin zone 


honeycomb lattice 


(direct space) (reciprocal space) 


Fig. 3. Graphene’s honeycomb lattice and the first Brillouin zone. 


The first Brillouin zone of the reciprocal lattice is shown 
in Fig. 3. The reciprocal space is a convenient concept for 
describing the reciprocal quantities (usually momentum) in 
the Fourier domain. It is defined by the reciprocal basis 


3, (8-55), 


b, = <*(x+ N35) and b, = 


which are related to direct vectors a the orthogonal 


relation b, = 210, . The reciprocal vectors are 
constructed from the direct space by the relations 
zZxa zZxa 
b, S= 2a 2 5 b, = 24 ——_ (1) 
Z-a xa, Z-a xa, 


The six points at the corners of the Brillouin zone are called 
the Dirac points. Only two of these points, represented as 


K = (Vss) and K’ = or, (v3 F). 


are inequivalent, and the rest are their periodic images 
connected by the reciprocal lattice vectors. The nearest 
neighboring atoms are represented by the vectors 


8, => (%+ v39). 8, =5(%- v3y) and 0, =— 


These vectors are used to derive the tight binding energy 
band structure in the next section. 








C. Energy band structure (tight binding approach) 


Tight binding is an efficient method for approximating the 
band structure, using superposition of a set of basis functions 
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closely related to the electronic orbitals of isolated atoms. In 
this model, it is assumed that the crystal potential only 
weakly perturbs the electronic orbitals of the isolated atoms. 
Therefore, the wave function of the electrons in the crystal 
structure can be approximated by the linear combination of 
the wave functions of the electrons in the isolated atoms. 


Since the overlap of the 2p, orbitals in graphene is small, 


this method provides a good description of the band structure 
of graphene. The wave function of the electron is then 


described as the linear combination of the 2p, orbitals over 


the entire crystal. It may be written as the superposition of 
the orbitals of the carbon atoms at the A and B sites, 
namely 


W=CiP, + CeQp (2) 


where @, and @, are formed, using the Bloch theorem, as 


the linear combination 


9, => exp(jk-r,)X(r-r,) (3) 
A 

0, = > exp(jk-r,)X (r-r,) (4) 
B 


where X (r) represents the wave function of the 2p, 


orbital. Due to the small overlap between the 2p, orbitals 


between neighboring atoms, the overlap integral is negligible 
in the tight binding development [17], and may, therefore, be 
dropped, 1.e. 


[Xr -r,)X (r -r;)dr = 0 (5) 


The energy states are found by solving the Schroedinger 
equation 


Hy = Ey (6) 


where H is the Hamiltonian of the crystal, which can be 
written as H=H,+AU, with HA, representing the 


Hamiltonian of isolated carbon atoms at their sites and AU 
representing the perturbation by the neighboring atoms. 
Replacing yw in (6) by (2), successively multiplying both 


sides by p; and P; and integrating the resulting equations 


over the entire lattice yields [17] 


jC, +H ypc, = ESC, 


(7) 
HC, + HypC, = ESC, 


where 


H, =| ¢\He,dr. H s = Hy = | of odr, 
Hy, =|9;H@,dr, S=|ppdr=| opdr ©) 


The energy eigenvalues are found by nullifying the 
determinant of (7), 


H -ES H 


AB ~() (9) 
A ,, Hm= ES 


The assumption (5) leads to S = N , where N is the total 
number of atoms in the crystal, and where the crystal 


symmetry implies H ,, = H,,. Therefore (9) becomes 





l / / 
E= yy (Aad HA NS E (10) 
with 
pad d = 
5 = PLT are) | (r-r,)HX(r-r,)dr 
(11) 
bi cot . š a 
U ATE —r,)|[X (r—r,) HX (r-r, )dr 
(12) 


Keeping only the nearest neighbour contributions in (11) and 


(12), H’, and H’; are approximated as 


Hig BorN Š exp(jk -a,) 


i=1...3 
= E, + yal cos(V3k,a)+ sco? ka SED) 


H a= > exp(—jk -6, ) 


i=1...3 


=o). op (ika) +2 cos Baa) an(-i4ka) 


(13) 
where Èp is the energy of the electron in the 2p, orbital of 
an isolated carbon atom [17], 7,, is the nearest neighbor 


hopping energy with an estimated value of 2.8 eV, and 7, » 


is the next nearest neighbor hopping energy with an 
estimated value of 0.1 eV, which is much smaller compared 


to Y¥4,. The energy band structure E(k) = E(k) SE AS 
thus given by 
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3 9 3 ; 
e(K+q) = OY py eS rgalal-[S 7 +—y,.a° sin (30, Ja Juri 
e(k)=+7 943+ f(k)+74,f(k) 2 4 2 


q 
3 3 0 = arctan| ~ j 
f(k) =2cos(V3k,a)+ 4cos Beg co a q q. (15) 
(14) 
[4], plotted in Fig. 4. Therefore, neglecting the terms with the factor 7,,, the 


band structure at the Dirac points is approximated by the 
linear symmetric Dirac cone equation 


e(q)= thv,q| (16) 


3 = 
where v, = pa '~10° m/s is the Fermi velocity, 


and q is the momentum measured with respect to the Dirac 


points. The striking feature of this energy dispersion function 
is that the velocity of the carriers is independent of their 
energy, similar to dispersion of a photon in free space. In this 
sense, the carriers in graphene are called massless Dirac 
fermions. Unlike massive particles, whose velocity is 
increased as the square root of the energy, electrons in 
graphene are, energy-wise, instantly accelerated to the Fermi 
velocity, which is the mechanism behind its high mobility. 
Another interesting property of the Dirac cone is the electron 
hole energy symmetry, which is the reason for graphene’s 
ambipolarity. The high symmetry of the band structure also 
leads to strong interband transition, making graphene a good 
candidate for light modulators [18-20]. Since graphene has a 
gapless band structure, the interband conductivity can be 
easily tuned for terahertz, infrared or optical frequencies by 
applying the appropriate amount of electrical doping, with 
applications in tunable graphene-based components such as 
optical modulators [18-20]. Another interesting property 
associated with the Dirac cones is the ability to provide 
unusual conductivities. For some specific parameters, the 
gapless linear band structure leads to conductivity with a 
dominant interband contribution. This condition results in a 
conductivity with a positive imaginary part, as opposed to 
the Drude conductivity, leading to exotic effects such as TE 
surface plasmons [5]. 





Fig. 4. Energy band structure of graphene for Y4g = 2.8 eV and D. Dirac Hamiltonian 


Nis ~ 0.1 ev. (a) Perspective plot of the conduction and valence bands, Writing (7) in the following form 


(b) projected color plot of the conduction band. The k. and k y axes have 


H' H' 
dimension of | / @ . The energy is in eV. aA an y = Ewy (17) 
H, H; 
BA BB 


The unique properties of graphene’s band structure is 
revealed by inspecting the Dirac points. Expanding (14) in 
Taylor series around the Dirac points K (or K’) for 


k= K +q for q << K gives [4] 


reveals that this equation satisfies the Schroedinger equation 
with the new 2-spinor wave function 
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C4 

Y= | (18) 
Cp 

and the effective Hamiltonian 
H'. #' 

Hel a (19) 
Al BA Al BB 


where the two components of the wave function are 
respectively connected to the A and B sub-lattices through 
(2). Expanding the effective Hamiltonian (19) around the 
Dirac point K results in 


0 q4,- jq, 


H, =i 
oe qd, + J4, 0 


=hv,o-q (20) 


satisfying H, w =e(q)y, where 6=(0,,0,) is a 


vector of Pauli matrices [4]. Equation (20) is the Dirac 


equation for massless spin- > particles (Weyl equation), 


operating on 2-spinor wave functions, with the velocity of 
light replaced by the Fermi velocity. In analogy with 2- 
spinors describing real spin, the 2-spinor wave function (18) 
describes what in literature is referred to as pseudo-spin and 
originates from the existence of two  sub-lattices in 
graphene [1]. The effective Hamiltonian around the K’ 
point is similarly expressed as 


H,. =hv,6"-q (21) 


II. GRAPHENE UNDER A STATIC MAGNETIC BIAS 


First, graphene’s effective Dirac Hamiltonian is used to 
describe the quantized energy levels of magnetically biased 
graphene. Then, based on transitions between these levels, 
the Kubo conductivity tensor is presented and shown to 
reduce to a Drude conductivity tensor under conditions that 
typically prevail in the devices to be described in Section V. 


A. Landau levels 


Under a magnetic field bias, the carriers are deflected by the 
Lorentz force and move in cyclotron orbits. The magnetic 
bias modifies the Hamiltonian of the crystal and limits the 
allowed cyclotron energies to a set of quantized Landau 
levels. In conventional parabolic energy dispersions, 
occurring in conventional two-dimensional electron gas 
(2DEGs), these levels are equidistant, with an energy 


difference ho,, where @. is the cyclotron frequency [21]. 


For relativistic linear band structures, the quantized allowed 
energies take a completely different form. Consider the Dirac 


Hamiltonian for graphene under a normal static magnetic 
field 8 . The states near the K point satisfy the Dirac 
equation corresponding to (20), 


hv ,6- QW, =EW, (22) 
or, explicitly, 


à pe 
Tee 
q, + jq, 0 Wp Wke 


Eliminating Y,, from the two equations results in 


hv, 


(Av, y (a, + ją, (a, - ją, )V xe = EY kp (24) 


The momentum operator is hq =—ihV — eA , where A is 
the magnetic vector potential. We will use the Landau gauge 


for the magnetic vector potential, i.e. A =(—By,0,0). 


Replacing the momentum operators hq. = —ih0 — eA. 
and hq, = -iho , reduces (24) to 


Ve ((-ind, + eBy) + eBh— no, Wes (x,y) = EY kp (x,y) 
(25) 


The wave function general solution to this equation is 


W kp (x,y) =e p(y) leading to 


Lea) (y+ Ms) -pea Joo) = (2) “re 


(26) 


which is the well-known equation of a quantum harmonic 


ik „x 


oscillator with frequency @ = leB | and mass m=1 [22]. 


The discrete set of the eigenvalues of (26) are given by 


2 
d na +ha =ho( n>) 
ZN 2 


which, upon solving for ¢, finally gives the following set of 
quantized Landau levels 


(27) 


E = EV 


n 


2h|\eB| n (28) 
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where n is an integer. Therefore, unlike in classical 2DEGs, 
the quantized Landau energy states in graphene are 
proportional to the square root of the magnetic field and 
spaced unequally. As the energies increase, the distance 


between the Landau levels decreases as €,,, —&, © ——. 

Vn 
Consequently, for high enough energies, these levels become 
very close to each other. Eventually, these discrete energy 
levels can be treated as a continuum, leading to the 
simplification of the conductivity models described in the 
next section. 


B. Conductivity model 


Interaction of the carriers in graphene with a tangential 
electric field is modelled through the conductivity given by 
the Kubo formula. The Kubo formula describing general 
many body problems is derived in the context of linear 
response of the quantum mechanical many body system to a 
small perturbation [23]-[25]. For an energy-independent 
scattering rate, the graphene’s Kubo conductivity is reduced 
to the following useful Lorentzian form, giving the diagonal 
and off-diagonal components of the conductivity as [26]- 
[29] 


o, (a, u,,B joe 
d 5 FCF 09 ae 
= (M,.-M,) -(@- 727) 
x ;-_ A —— 
M Ma MAM, 
(M +M,) -(@- jT y h? 


n+l 


(29) 


27 
o, (©, H., BoT) =e 


o0 


xy (M,)- f (Ma) SA Malr (-M,)} 


n=0 


A? l 
Jit 
MM )(M,.—-M,) -(@-j2r) nh 
er aay. 
le eo oi ee oe ne 
M,M mi (M, +M,) -(@- j2) h 
(30) 


where 4, is the chemical potential, I =1/2r with T 


-1 
; : — u, )/kgT ; 
represents the scattering time, Ja = (I+e" HeMkp 1S 


the Fermi-Dirac distribution, and 


M = JA? + hiv; leB,|n is the m-th Landau level 


energy. It has been speculated that under very high magnetic 
fields, a gap might open up in the energy band structure [26]- 
[28], called an exitonic gap described by the parameter A. 
For typical values of the magnetic field considered in this 
paper, A is set to zero. Equations (29) and (30) involve 
energy transitions within the conduction, or valence band 
called the intraband transitions, and energy transitions 
between valence and conduction bands called interband 
transitions. For a pedagogical description of these energy 
transitions the reader is referred to [30]. In the low magnetic 


field limit B,—>O , (29) and (30) have a poor 


convergence. However, since the Landau levels become very 
dense these series converge into integrals more appropriate 
for numerical purposes (see [29]). 


In the limit of no interband transition, 4, >> 2h@, weak 
Landau quantization (dense Landau levels), 44, >> ho, , 


and single carrier type conduction, u, >> k T , the Kubo 


conductivity reduces to the classical Drude form 


1+ j@T 
o,(a,B, )=o, ———— (31) 
(mB) =o (w,r) +(1+ jor) 
OT 
o, (a, B, )=o, ———_; (32) 
(2, By) =o (ar) +(1+ jar) 
with the DC conductivity 
o 2er T In 2. cosh (33) 
"gh? ” 2k, 


The conductivity tensor is thus asymmetric. An important 
property of an anti-symmetric conductivity tensor is 
discrimination between right- and left-handed circularly 
polarized excitations. An anti-symmetric conductivity tensor, 


o=0,(xxt+yy)+o,(yx—xy), has two circularly 
x+y 


eigenvalues O} Ł JO . The response of such conductivity 


polarized eigenvectors corresponding to the 


to a circularly polarized electric field E}, = E (xF jy) 


can be expressed as 
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F — 
‘Kop — 


Qıl 


+ JY 
J 


o,)(x¥ jy) (34) 


E,o-(x 
L (o; £ 
(o, + jo,)E 


A H Ne NN 


P 


Equation (34) reveals that the conductivities seen by 
right/left-handed circularly polarized waves are the scalar 


conductivity O} + Jo,. As a result right- and left-handed 


circularly polarized waves see different material parameters 
and behave differently. This process is the cause of 
nonreciprocity in the magnetoplasmons to be discussed in the 
next sections. 


IV. GRAPHENE MAGNETOPLASMONICS 


Surface plasmons are collective oscillations of charges at 
the interface of metal and dielectric media or a 2DEG, as 
represented in Fig. 5. The electromagnetic field created in 
the process exponentially decays from the surface and is thus 
tightly confined to it. The perturbations produced by the 
collective oscillation of charges create an electric field in the 
direction of propagation, while the magnetic field is 
transverse to the k vector. These types of surface waves 
which propagate in conventional 2DEGs and in graphene are 
referred to as longitudinal or transverse magnetic (TM) 
surface plasmons. TM surface plasmons can propagate in 
materials whose conductivity take a negative imaginary part 


. lot .. . 
(assuming thee” time convention). 







positive charges 


negative charges 


Fig. 5. TM surface plasmons, consisting of longitudinal electric fields due 
to the longitudinal displacement of charges at a metal-dielectric interface or 
on a 2DEG. 


TM surface plasmons in 2DEGs (including graphene) 


exhibiting a Drude conductivity response have a @ © Vk 
dispersion and are, therefore, characterized by a wavelength 
that becomes highly compressed and an energy that becomes 
tightly confined to the surface at high frequencies. Such 
surface plasmons exhibit lower loss at frequencies where the 
imaginary part of the conductivity increases to a point where 
it becomes comparable to its real part. In the Drude regime 


this condition occurs at a frequency corresponding to the 
scattering rate. For graphene, this frequency window starts at 
terahertz frequencies and, due to lower scattering, graphene 
plasmons have lower loss compared to other 2DEGs. 


Unlike other 2DEGs, graphene’s conductivity can take a 
positive imaginary part under specific conditions. Such 
conditions occur at frequencies where the interband 
conductivity becomes dominant [29]. In this case, graphene 
supports transverse electric (TE) surface plasmons (referred 
to as transverse or TE or H surface plasmons) [5]. TE surface 
plasmons are less dispersive and less tightly confined 
compared to their TM counterparts, and have a phase 
velocity close to the speed of light in the surrounding 
medium. However, they exhibit lower loss compared to TM 
modes. 


Magnetoplasmons are surface charge oscillations in the 
presence of a static magnetic field. Under the influence of 
the magnetic field, carriers move in cyclotron orbits and, as a 
consequence, all the components of the electric and magnetic 
fields are present. Therefore, magnetoplasmons in infinite 
graphene are hybrid (TM-TE) modes. The dispersion relation 
of magnetoplasmons in infinite graphene is obtained by 
applying electromagnetic boundary conditions on the 
graphene layer to a surface wave. Assuming an infinite 
graphene sheet in the xz plane and a surface wave 
propagating in the z direction, with an electric field 


— jk, z— ‘ : : 
E = Ee kee arp] , the dispersion relation reads [8] 


2 72 2 
k; =k +a; 


2 
L482 +s? t,/(14+52+5?) —452 O9 
. d o d o d 


a 
Y= 
ko : 28); 


where S$, =79,0,/2, 8, =O, /2. 


In a graphene strip with a finite width, plasmons or 
magnetoplasmons exhibit transverse resonances across the 
strip. The corresponding modes are called the 2D bulk modes 
of the graphene strip. A magnetically unbiased graphene 
strip supports an infinite number of bulk modes and two 
degenerate edge modes, each propagating on a different 
edge. The dispersion of these modes are shown in Fig. 6 
with solid red and blue curves. The dispersion of 
magnetoplasmons in an infinite graphene sheet with the same 
material parameters is plotted in a dotted line for 
comparison. 
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Fig. 6. Dispersion curves for magnetoplasmons in a graphene strip with 
parameters: width w = 100 um, relaxation time t= 0.1 ps, carrier density 
i= 10° cm” and B, = 1 T. Edge modes are plotted in red and bulk 


modes in blue. The dotted curve shows the dispersion for an infinite 
graphene sheet with the same parameters. 


The edge modes have exponentially decay away from the 
edge, and therefore have an imaginary transverse 
wavenumber. Considering the strip dispersion relation, 


ke = ki -k -k = ki +a? +a , and taking into 


account that magnetoplasmons in infinite graphene have 
exponential decay only along y, the strip modes are expected 


to have a larger k, , which explains why they lie beneath the 
infinite graphene line in Fig. 6. The strip bulk modes have a 


sinusoidal variation in the x direction with a real k. and 
therefore, have a smaller k, compared to infinite graphene 


with k, =(), lying above the infinite graphene dispersion 
curve as expected. 


A closer look at the electromagnetic fields of the edge 
modes in the plane of the graphene strip reveals an 
interesting feature of their polarization. The electric field 
inside the plane of the strip for the left and right edges are 
shown in Fig. 7. The vector field plots show the two edge 
modes are circularly polarized with opposite (left- and right- 
handed) polarizations. In the presence of a static magnetic 
bias, the right-handed and left-handed circularly polarized 
waves see different conductivities according to (34). 
Therefore, application of a magnetic bias lifts the degeneracy 
of the edge modes and confers in them different dispersions. 
By a similar token, it may be inferred that edge modes of the 
two edges show asymmetric dispersions for the opposite 
directions of propagation. This phenomenon is exploited in 
the next sections in the design of nonreciprocal 
magnetoplasmonic devices. 








(b) 


Fig. 7. Electric field on the graphene strip for the edge modes 
propagating on the right and left edges. (a) Point A on the right edge sees a 
counter clockwise rotating electric field as the wave (mode 2+) propagates 
along the graphene strip. (b) Point B on the left edge sees a clockwise 
rotating electric field as the wave (mode 1+) propagates along the graphene 
Strip. 


V. NONRECIPROCAL GRAPHENE MAGNETOPLASMONIC 
DEVICES 


A. Nonreciprocal phase shifter 


Under a static magnetic bias, the edge magnetoplasmons 
in a graphene strip exhibit different dispersions for opposite 
directions of propagation, as seen from the asymmetry of the 
thick red line (left edge mode) in Fig. 8, and similarly for the 
right edge mode. 
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Fig. 8. Dispersion curves for magnetoplasmons in a graphene strip with 
parameters: width w = 100 um, relaxation time t= 0.1 ps, carrier density 
n, = 10° cm” and B, = 1 T. Edge modes are plotted in red and bulk 
modes in blue. The dashed curve shows the dispersion for an infinite 


graphene sheet with the same parameters. The gray area corresponds to the 
light cone. 


One of the consequences of such a nonreciprocity is that the 
edge modes produce different phase shifts in the forward and 
backward directions. Assume that the right edge mode is 
shorted and the left edge mode is excited, as illustrated in 
Fig. 9. Due to the dispersion asymmetry, the excited edge 
mode will produce different phase shifts in opposite 
directions, thus realizing a nonreciprocal phase shifter. The 
phase shifts, plotted in Fig. 10 for two different carrier 
densities, show a significant phase difference between the 
opposite directions. 





Fig. 9. Graphene-based nonreciprocal phase shifter. The right edge mode 
is shorted by a PEC. The left edge mode produces different phase shifts in 
forward and backward directions. 
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Fig. 10. Phase difference and loss produced by a nonreciprocal phase 
shifter for carrier concentrations n = 10'°° cm™~ and n, = 5x10? cm” 


B. Nonreciprocal coupler 


The nonreciprocity of edge magnetoplasmons can also be 
used for the realization of exotic nonreciprocal components, 
such as nonreciprocal plasmonic couplers. The coupler 
structure is depicted in Fig. 11 [9]. It consists of two in-plane 
edge-coupled magnetically biased graphene strips. The 
coupler is designed such that for the forward direction, the 
edge mode excited at port 1 couples to port 4. In the 
backward direction (excitation at port 2) however, coupling 
does not occur. 


(a) (b) 


Fig. 11. Non reciprocal plasmonic coupler, consisting of two parallel 
graphene plasmonic waveguides. Both waveguides are biased with a 
magnetostatic field perpendicular to their plane. (a) Feeding at port 1. 
(b) Feeding at port 2. 


To design such a coupler, the physical parameters of the 
graphene strips should be tuned in such a way that in the 
forward direction, the edge modes of the adjacent edges of 
the strips are phase matched and in the opposite direction 
their phase constants must be unequal. To achieve this 
condition, the phase constants are adjusted by tuning the 
carrier density in the strips. The phase constants of the 
magnetoplasmon modes of the two isolated graphene strips 
are shown in Fig. 12. The modes of the right strip are shown 
in solid curves and the modes of the left strip in dashed 
curves. Red colors represent the edge modes and blue colors 
the bulk modes. 
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Fig. 12. Dispersion curves for edge (red) and bulk (blue) modes of two 
isolated graphene strips with different parameters. The solid curves show the 
slow-wave factor and loss for a graphene strip with parameters, width w = 
100 um, relaxation time t= 0.1 ps, carrier density n, =10! cm” and 


B= 1 T. The dashed curves show the slow-wave factor and loss for a 


graphene strip with parameters, width w = 100 um, relaxation time t= 
0.1 ps, carrier density —8x10' cm” and B, = 1 T. Phase matched 


regions are emphasized by ellipses. 


The electric field of the different edge modes of isolated 
strips are plotted in Fig. 13, showing different possible 
scenarios for the interaction of the edge modes propagating 
on the left and right strips. In the operation zone marked by 
the large ellipse in the forward direction, the edge mode 


propagating on the right edge of the left strip has a k, vector 


matched to the phase constant of the edge mode propagating 
on the left edge of the right strip. This scenario is highlighted 
in Fig. 13 by the red rectangle. 


phase 
matched mismatched 


phase 





left strip 
forward (+z) 


right strip 
backward (—z) 


Fig. 13. Electric field magnitude for the edge modes of the graphene strips 
of Fig. 11 showing different possible scenarios when they are placed side by 
side. The right strip has parameters with parameters: width w = 100 um, 
relaxation time t= 0.1 ps, carrier density n, = 10° cm and B ITI 


The left strip has parameters with parameters: width w = 100 um, relaxation 
time t= 0.1 ps, carrier density n = 8x10 cm” and b~ lI 


Thus, when the two strips are close to each other, the 
mode excited at port 1 couples to port 4. In the backward 
direction, the situation is different. The modes shown in the 
highlighted region by ellipses are phase matched. However, 
they are propagating on the opposite edges, as indicated in 
Fig. 13 by the green rectangle. The modes propagating on the 
inner edges, shown in Fig. 13 by the blue rectangle, are 
phase mismatched, as seen from their dispersions in Fig. 12 
indicated by small circles. 


Figure 14 shows the simulated electric fields of the 
coupler in the forward and backward directions. In the 
forward direction, the edge modes couple and form 
symmetric and anti-symmetric modes with different phase 
constants, thus realizing a forward coupler. In the backward 
direction, however, the edge modes do not couple. 


. 4% Ea n 
(a) forward (+z) (b) backward (—z) 


Fig. 14. Electric field magnitude for the edge modes of the coupler of 
Fig. 11 propagating on the near edges of the strips at the frequency f = 6 
THz. (a) Forward direction. The edge modes couple and give two coupled 
symmetrical and anti-symmetrical modes. (b) backward direction. The 
modes do not couple. 


The dispersion curves for the modes of the coupler are 
plotted in Fig. 15. The symmetrical and anti-symmetrical 
dispersion curves are plotted in black. The symmetrical and 
anti-symmetrical modes have different phase constants, and 
therefore, the edge mode excited at port 1 in the forward 
direction will couple to port 4. 





2-11 2 
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Fig. 15. Dispersion curves of the edge and bulk modes of the graphene 
plasmonic coupler of Fig. 11 for the forward and backward propagation. 
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The vector plot of the transverse electric field for the 
symmetric and anti-symmetric modes are shown in Fig. 16. 
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Fig. 16. Transverse electric field vector plot for the coupled edge modes. 
(a) anti-symmetric and (b) symmetric. 


C. Plasmonic isolators 
1. Chemically doped isolator 


Another interesting nonreciprocal plasmonic mode is the 
mode propagating along a graphene p-n junction. A 
chemically doped graphene p-n junction structure is depicted 
in Fig. 17, consisting of two chemically p-doped and n- 
doped in-plane graphene strips separated by a nano-gap [12]. 
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Fig. 17. Chemically doped graphene p-n junction. 


The dispersion of the magnetically unbiased structure is 
shown in Fig. 18. This structure supports two edge modes, a 
p-n junction mode, shown in red, and an infinite number of 
bulk modes propagating inside the strips. 





1 15 i 2.5 4 


—Im (k) ko) 


Fig. 18. Slow wave factor and loss for the structure in Fig. 17 with no 
magnetic bias. w= 100 um, t= 0.1 ps, n= p=10" cm~ , T= 300 K and 
B~ 0 T. The p-n junction mode is plotted in red. 


Inspection of the electromagnetic fields in the plane of the 
strips reveals interesting facts about the p-n junction mode. 
The electric field of this mode is shown in Fig. 19 for the 
unbiased structure. In the forward direction, point L on the 
left strip sees a right hand circularly polarized (RHCP) wave 
under propagation, and point R on the right strip sees a left 
hand circularly polarized (LHCP) wave. As the magnetic 
bias is switched on, the conductivity becomes tensorial and 
RHCP and LHCP waves see equivalent conductivities of 
O; + jO, and O} — JO, , respectively. However, the p and 
n doped strips will exhibit off-diagonal conductivities with 
opposite signs O,, =—O,,. 


of the junction 


Therefore, in the forward 


direction, both sides see the same 


conductivity, namely O} + JO. Following a similar logic, 
in the backward direction, the waves on both sides see the 
same conductivity, namely O} — JO. Therefore, in the 


presence of a magnetic bias, the p-n junction mode behaves 
differently in the forward and backward directions. In the 
forward direction, it becomes concentrated at the junction as 
the magnetic field is switched on. In the backward direction 
it splits away from the junction. The dispersions of the 
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modes and the electric field patterns of the junction mode are 
plotted in Fig. 20. 
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Fig. 19. Electric field pattern in the plane of graphene for the structure of 
Fig. 17 with B ,=0 T. Points R on the right strip sees a clock-wise rotating 


electric field in the forward direction and a counter clock-wise rotating 
electric field in the backward direction as the wave propagates. Point L on 
the left strip sees an oppositely rotating electric field to point R in each 
direction. 
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Fig. 20. Slow-wave factor and loss for the magnetoplasmonic isolator in 
Fig. 17;. w= 100 mm, t= 0.1 ps, n= p=10" cm~ , T = 300 K and B, = 


1 T. The p-n junction mode is represented in red. 


The evolution of the p-n junction electric field as the 
magnetic field is increased is shown in Fig. 21. This 
nonreciprocal phenomenon can be exploited for the 
realization of plasmonic isolators. If the magnetically biased 
structure is excited at the center with the appropriate field 
distribution, in the forward direction, it will couple to the p-n 
junction mode. However, in the backward direction, it will 





be reflected back as there is no mode with energy at the 
center. 
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Fig. 21. Evolution of the electric field pattern of the mode propagating at 
the p-n junction (mode 1) for the forward and backward propagations as the 
magnetic field is increased. 


2. Electrically doped isolator 


One of the benefits of graphene is the possibility of tuning 
the carrier density through electric biasing. An alternative 
method for producing p-n junctions in graphene is through 
the application of a tangential electric field. 
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Fig. 22. Electrically doped graphene p-n junction. 


The electric field changes the balance of the carriers on the 
strip and creates p- and n- doped regions on the right and left 
sides of the strip. The density of the carriers is found by 
solving the 2D Poisson equation. The unknown charge 


density p(x) on the strip is determined by solving the 


integral equation 


12 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


w/2 


| (x, y')G(x,y;x', y')dx' -Ex =0 
—w/2 
—w/2<x<w/2, 


(36) 


y=0, y'=0 


where the 2D Green function solution to the Poisson 
equation is given by 
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Fig. 23. Net carrier density and electric potential for a graphene 
strip doped with an electric field for w= 50 um and £ =10° V/m. 


The resulting carrier density is shown in Fig. 23 for a 
graphene strip with parameters w= 50 um, £ =108 V/m, T= 
300 K. Note that although the net carrier density is zero near 
the center, the thermally excited carriers produce a 
significant conductivity at the junction. To find this 
conductivity, one first computes the chemical potential 
across the graphene strip by solving the integral equation 


0 


i. =f (2,u,)N(e)de-[[1-f,(e,u,)|N(e)de 


—0oo0 





2\é 
n(e)= ŻE G8) 
mh vi 
where f, is the Fermi-Dirac distribution and M, „ is the net 


carrier density, plotted in Fig. 23. The resulting density is 
plotted in Fig. 24. Note that in the Drude regime the 
conductivity is proportional to n+ p, where n and p are 


the density of electrons and holes respectively. 
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Fig. 24. Thermally excited electron and hole densities at room 


temperature for a graphene strip doped with an electric field; w = 50 um, 
E, =108 V/m, T= 300 K. 


Finally the non-uniform conductivity across the strip is found 
by substituting the chemical potential inside the Kubo 
formula. This magnetically unbiased structure supports a p-n 
junction mode with maximum energy at the junction, an 
infinite number of bulk modes, shown in Fig. 25, and two 
edge modes (not shown here). 


J (THz) 





Fig. 25. Slow-wave factor and loss for a graphene strip biased by an 
electric field; w = 50 um, Le 10° V/m, T = 300 K. The p-n junction 


mode is represented in red. The insets show the electric field pattern. 


In the presence of a magnetic bias, the p-n junction mode 
sees different conductivities in the forward and backward p 
propagation directions, as understood by following a similar 
argument as for the chemically doped case, and therefore 
becomes nonreciprocal. In the forward direction, this mode is 
concentrated at the center and, in the backward direction, it 
splits away from the center, as shown in Fig. 26. This 
phenomenon can be exploited for the realization of tunable 
plasmonic isolators. 


13 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


f (THz) 





-50 -40 -30 -20 -10 10 20 30 40 +0 


Re(k-/ ko) 


Fig. 26. Dispersion curves for a magnetically biased graphene strip 
biased by an electric field; w = 50 um, E, =108 V/m, T = 300 K, 
B,=0.1T. The p-n junction mode is represented in red. The grey area 


represents the light cone. 


VI. CONCLUSIONS 


First, the fundamentals of graphene, including carbon 
bonds, crystallographic characteristics, energy band 
structure, Dirac Hamiltonian, Landau quantization and 
conductivity model, was introduced in a self-consistent and 
pedagogical fashion. Next, graphene magnetoplasmonic 
principles were overviewed. Finally, the nonreciprocal 
properties of magnetoplasmonic modes were described and a 
number of related device applications, such as isolators and 
couplers, were presented. 
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